Electrothermal feedback in superconducting nanowire single-photon detectors 
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We investigate the role of electrothermal feedback in the operation of superconducting nanowire 
single-photon detectors (SNSPDs). It is found that the desired mode of operation for SNSPDs 
is only achieved if this feedback is unstable, which happens naturally through the slow electrical 
response associated with their relatively large kinetic inductance. If this response is sped up in an 
effort to increase the device count rate, the electrothermal feedback becomes stable and results in 
an effect known as latching, where the device is locked in a resistive state and can no longer detect 
photons. We present a set of experiments which elucidate this effect, and a simple model which 
quantitatively explains the results. 
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Superconducting nanowire single-photon detectors 
(SNSPDs) p], 0, y, El have a combination of attributes 
found in no other photon counter, including high speed, 
high detection efficiency over a wide range of wave- 
lengths, and low dark counts. Of particular importance 
is their high single-photon timing resolution of ~30 ps 
[3], which permits extremely high data rates in commu- 
nications applications [5|, [6[. Full use of this electrical 
bandwidth is limited, however, by the fact that the max- 
imum count rates of these devices are much smaller (a 
few hundred MHz for 10 pm 2 active area, and decreas- 
ing as the area is increased Q), limited by their large 
kinetic inductance and the input impedance of the read- 
out circuit 0, 0- To increase the count rate, therefore, 
one must either reduce the kinetic inductance (by us- 
ing a smaller active area or through the use of different 
materials or substrates) or increase the load impedance 
Q. However, either one of these approaches causes the 
wire to "latch" into a stable resistive state where it no 
longer detects photons [8]. This effect arises when nega- 
tive electrothermal feedback, which in normal operation 
allows the device to reset itself, is made fast enough that 
it becomes stable. We present experiments which probe 
the stability of this feedback, and we develop a model 
which quantitatively explains our observations. 

The operation of an SNSPD is illustrated in Fig. Ufa). 
A nanowire (typically ~100nm wide, 5nm thick) is biased 
with a DC current io near the critical current I c . When 
a photon is absorbed, a short (<100 nm long) normal do- 
main is nucleated, giving the wire a resistance R n (t); this 
results in Joule heating which causes the normal domain 
(and consequently, R n ) to expand in time exponentially. 
The expansion is counteracted by negative electrothermal 
feedback associated with the load Rl (typically a 50O 
transmission line), which forms a current divider with the 
nanowire and tends to divert current out of it, thereby 
reducing the heating (and producing a voltage across the 
load). In a correctly functioning device, this feedback is 
sufficiently slow so that it is unstable: that is, the induc- 
tive time constant is long enough so that before apprecia- 
ble current is diverted out of the load, Joule heating has 



already caused the normal domain to grow large enough 
so that R n > R L [ll]|. The current in the device then 
drops nearly to zero, allowing it to quickly cool down 
and return to the superconducting state, after which the 
current recovers with a time constant r e = L/Rl [2|. If 
one attempts to shorten r e , at some point the negative 
feedback becomes fast enough to counterbalance the fast 
Joule heating before it runs away, resulting in a stable 
resistive domain, known as a self- heating hotspot 0. 

In a standard treatment of these hotspots @ , solutions 
to a one-dimensional heat equation are found in which 
a normal-superconducting (NS) boundary propagates at 
constant velocity v ns , for fixed device current Id E3] • 
This results in a solution of the form: 
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where = V 'A cs nh/c is a characteristic velocity (A cs is 
the wire's cross-sectional area, k is its thermal conduc- 
tivity, and c and h are the heat capacity and heat trans- 
fer coefficient to the substrate, per unit length, respec- 
tively), I c is the critical current, and a = p n I 2 /h(T c — To) 
is known as the Stekly parameter (p n is the normal 
resistance of the wire per unit length, To and T c are 
the substrate and superconducting critical temperatures) 
which characterizes the ratio of joule heating to con- 
duction cooling in the normal state [9]. Equation [T] 
is valid when To <C T c , and the approximate equal- 
ity holds for sm all devia tions SI = I — I ss , with 7 = 
(T c - T )(c/p n )^/h/^A cs and I 2 S = 2h(T c - T )/p n . The 
physical meaning of eq. [T] is clear: the NS boundary 
is stationary only if the local power density is equal to 
a fixed value; if it is greater, the hotspot will expand 
(v ns > 0), if less it will contract (v ns < 0). 

We can use eq. Q] to describe the electrothermal cir- 
cuit in fig. [2(a), by combining it with: dR n /dt = 2p n v ns 
and IdRn + Ldld/dt = Rl(Io — Id)- Assuming that the 
small-signal stability of the DC hotspot solution (Id = 
I SS ,R SS = Rl(Iq/Iss — 1)) determines whether it will 



latch, we linearize about this solution to obtain a second- 
order system, with damping coefficient £ = jf- yrthjr~e^ 
where r t h = Rl/^PuVo is a thermal time constant. This 
can be re-expressed as ( = \yJrth,tot/T e ,tou witn T e,tot = 
L/Rtot and r th ,tot = Rtot/2p n vo, where R tot = R L + R SS . 
If the damping is less than some critical, minimum value 
(latch, the feedback cannot stabilize the hostpot during 
the initial photoresponse, as described above, and the 
device operates normally. However, since the steady- 
state solution gives R ss oc io (as io is increased, a 
larger hotspot is necessary for Id = I ss ) the hotspot be- 
comes more stable as io is increased, until eventually 
((Iiatch) — (latch and the device latches. For a correctly 
functioning device, Iiatch > ic, so that latching does not 
affect its operation. However, if r e is decreased, iiatch 
decreases, and eventually it becomes less than I c . This 
prevents the device from being biased near I c |17| , result- 
ing in a drastic reduction in performance [12, LL3[ • 

Devices used in this work were fabricated from ~5 
nm thick NbN films, deposited on R-plane sapphire sub- 
strates in a UHV DC magnetron sputtering system (base 
pressure < 10 -10 mbar). Film deposition was performed 
at a wafer temperature of ~800 C and a pressure of 
~ 10 -8 mbar jl5j |. Aligned photolithography and liftoff 
were used to pattern ~100 nm thick Ti films for on-chip 
resistors [8], and Ti:Au contact pads. Patterning of the 
NbN was then performed with e-beam lithography using 
HSQ resist [3[. Devices were tested in a cryogenic probing 
station at temperatures of 1.8-12K using the techniques 
described in Refs. Q. 

Figure [2(c)- (f) show data for a set of (3/imx3.3/im ac- 
tive area) devices having various resistors Rs in series 
with the 50ft readout line [8[ [fig. [2(b)] so that Rl = 
50Q + Rs- Panels (c) and (d) show averaged pulse shapes 
for devices with Rs — and Rs — 250^, respectively. 
Clearly, the reset time can be reduced; however, this 
comes at a price. Panels (e) and (f) show, for devices with 
different Rs, the current /switch = min(i c , Iiatch) above 
which each device no longer detects photons, and the 
measured detection efficiency (DE) at Io = 0.975I sw u c h 
[13 • The data are plotted vs. the speedup of the reset 
time (Rs + 50^) /50O, and show that as this speedup is 
increased, I switch decreases (due to reduction of Iiatch), 
resulting in a significantly reduced DE. 

To investigate the latched state, we fabricated devices 
designed to probe the stability of self-heating hotspots 
as a function of io, and Rl- Each device consisted of 
three sections in series, as shown in Fig. [2] (a): a 3 pm- 
long, 100 nm-wide nanowire where the hotspot was nu- 
cleated [11[; a wider (200 nm) meandered section acting 
as an inductance; and a series of nine contact pads in- 
terspersed with Ti-film resistors. Also shown are the two 
electrical probes, which result in the circuit of Fig. [2](b): 
a high-impedance (R p — 20k^) 3-point measurement of 
the nanowire resistance. We varied Rl by touching the 
probes down to different pads along the line, and L by 
testing different devices (with different series inductors). 
We tested 66 devices on three chips, and selected from 
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FIG. 1: Figure [2 (color online) Speedup and latching of 
nanowire detectors with increased load impedance, (a) elec- 
trical model of detector operation. A hotspot is nucleated by 
absorption of a photon, producing a resistance R n in series 
with the wire's kinetic inductance L. (b) circuit describing 
experimental configuration for increasing load resistance us- 
ing a series resistor Rs. (c) and (d) averaged pulse shapes 
from 3/xmx3.3//m detectors with L ~ 50 nH for Rs — and 
Rs = 25011; red solid lines are predictions with no free pa- 
rameters, (e) switching current I SW itch = min(i c , Iiatch) vs. 
speedup 50Q/(Rs + 50Q). For small Rs, /switch = Ic but as 
Rs is increased, iiatch decreases, becoming less than I c . (f) 
DE at Io = 0.975/ sw itch vs. speedup (open squares). Also 
shown (crosses) are the expected DEs assuming that latching 
affects DE simply by limiting Io (obtained from DE vs. Io at 
Rs = 0). 



these only unconstricted [12] wires with nearly identi- 
cal linewidths (the observed I c of devices used here were 
within ~10% of each other - typically 22-24 /iA), span- 
ning the range: R L = 20 - 100012 and L = 6 - 600 nH. 

For each L and Rl, we acquired a DC I-V curve like 
those shown in Figs. [2](c),(d), sweeping Io downward 
starting from high values where the hotspot was stable 
[141 ]. All curves exhibit a step in voltage which can be 
identified as iiatch- For smaller Rl, where iiatch is large, 
I c can be identified as the point where an onset of re- 
sistance occurs common to several curves; this onset is 
more gradual with Io rather than sudden as at iiatch- 
These features can be understood by examining figs.[2](e) 
and (f), which show Id and R n inferred from the data of 
fig. [2](d). Above the discontinuous jump in current at 
^latch, Id is fixed at I ss (independent of io and Rl) in- 
dicating the latched state. For smaller r e , iiatch < Ic, so 
only iiatch is observed [l8[ . When r e is large enough that 
^latch > i c , an intermediate region appears where the re- 
sistance increases continuously with io; this arises from 
relaxation oscillations 

SGI, 

as indicated in the figure: 
the device cannot superconduct when io > i c , but nei- 
ther can a stable hotspot be formed when io < iiatch, so 
instead current oscillates back and forth between the de- 
vice and the load, producing a periodic pulse train with 
a frequency that increases as io is increased [17]. The av- 
erage resistance increases with this frequency, producing 
the observed continuous decrease in i^. 

Figure [3] shows the measured Iiatch as a function of Rl 
and L, which can be thought of as defining the boundary 
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FIG. 2: Figure [2] (color online) Hotspot stability mea- 
surements, (a) device schematic; two ground-signal-ground 
probes simultaneously perform a high-impedance 3-point 
measurement of the hotspot resistance, with Rl determined 
by probe position, (b) equivalent electrical circuit, (c) and 

(d) example V-I curves, with L = 60 nH and L = 605 nH, 
respectively; (e) and (f) inferred Id and R u /Rl for the data 
shown in (d). In (e) the relaxation oscillations in the region 
I c < Io < ^latch are shown schematically. Dashed lines show 

(e) I d = Io and (f) R ss = R L (I /I 8 s ~ 1). 



between stable and unstable hotspots. Our simple model 
described above predicts: r e /r t h cx (/i a tch/%s) 2 > a line of 
slope 1 in the figure (indicated by the dashed line). The 
data do approach this line, though only in the r e ^> r t h 
limit. This is consistent with the assumption of constant 
(or slowly varying) Id under which eq. [T] was derived. 
As r e /r t h is decreased, the data trend downward, away 
from this line, and I\&tch/Iss becomes almost independent 
of r e /rth] this implies a minimum I\ a tch/Iss, or equiva- 
lently, a minimum R ss /Rl, below which the hotpost is 
always unstable. This is shown in the inset: the mea- 
sured minimum stable R ss is always greater than ~ Rl. 

The crossover to this behavior can be explained in 
terms of a timescale r a over which the temperature profile 
of the hotspot stabilizes into the steady-state form which 
yields eq. [H For power density variations occuring faster 
than this, the NS boundaries do not have time to start 
moving, resulting instead in a temperature deviation AT. 
Since the NS boundary occurs at T w T c , where p n is 
temperature-dependent (defined by dp/dT = (3 > 0), this 
changes R ni giving a second, parallel electrothermal feed- 
back path which dominates for frequencies uo ^> r" 1 . We 
can describe this by replacing eq. [T]with: 
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FIG. 3: Figure [3] (color online) Summary of hotspot stabil- 
ity results. Data are shown from 3 different chips (indicated 
by different colors). Circles, squares, and triangles are data 
for L =6-12, 15-60, and 120-600 nH, respectively. In the 
r e ^> r a limit (where NS domain-wall motion dominates the 
electrothermal feedback), the data approach the dashed line, 
which is the prediction based on eq. [1] discussed in the text. 
For r e <C T t h the NS domain walls are effectively fixed and 
the temperature feedback dominates. In this regime the feed- 
back is always unstable when R ss ~> Rl (or equivalently, 
Io ~< 2I 38 ), as shown in the inset. The solid curves are 
obtained from eq. [8] assuming a phase margin of 30 degrees; 
each curve corresponds to a fixed L in the set (6,15,30,60,600) 
nH and spans the range of Rl in the data. The dotted lines 
extend these predictions over a wider range of Rl- 



teristic time over which dl/dt = 2v ns adapts to changes in 
power density: for slow timescales dt ^> r a , r a d 2 l/dt 2 <C 
dl/dt and eq. [2] reduces to eq. [T](with p = p n ). For faster 
timescales, r a d 2 l/dt 2 becomes appreciable, and acts as a 
source term for temperature deviations in eq. [3j When 
dt <C r a , r a d 2 l/ dt 2 ^> dl/dt and eqs. [2] and [3] can be com- 
bined to give: cdAT/dt w I 2 p- I 2 s p n - hAT 0. In this 
limit, if Rl ~ R ss the bias circuit including Rl begins 
to look like a current source, which then results in posi- 
tive feedback: a current change produces a temperature 
and resistance change of the same sign. Therefore, the 
hotspot is always unstable when R ss < Rl- 

Expressing eqs. [2] and [3] in dimensionless units (i = 
J d /Jo, r = R n /RL, A ee l(3T c /R L , = T/T c ), and ex- 
panding to first order in small deviations from steady 
state, we obtain: 
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Sr = ri(io-l)S9 + ri- 1 SX (5) 
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Here, I is the hotspot length, p(AT) is the resistance per 
unit length, and R n = p{AT)l. In eq. El r a is the charac- 
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Here, the prime denotes differentiation with respect to 
t/r e , i h/Iaa, 6 = (T c - T )/T c , rj pT c /p n char- 
acterizes the resistive transition slope, and r c = c/h 
is a cooling time constant. When r e ^> rt^,r a , the 
system reduces to: 8i" + io5i f — 4:T e /rth ~ 0, which 
has damping coefficient £ = io(4\/ r e/ r t/i) _1 ? as above. 
In the opposite limit, where r e <C r^,r a , we obtain: 
<K" + zo<K' + (2r]Qr e I V c ) (io — 2) « 0. In agreement with 
our argument above, the oscillation frequency becomes 
negative for R ss < R L (I < 21 3S ). 

We characterize the stability of the system of eqs. 0]- 
using its "open loop" gain A \\ we assume a small 
oscillatory perturbation by replacing Sr in eq. [H with 
Are ja;t , and responses (Si^O^SX^S^e^ 1 . Solving for 
A \ = Jr/Ar, we obtain: 
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The stability of the system can then be quantified by 
the phase margin: 7r+arg[^4 /(a;o)], where c^o is the unity 
gain (\A Q i\ = 1) frequency. In the extreme case, when the 
phase margin is zero (arg[A ^(o;o)] = — tt), the feedback is 
positive. The solid lines in Fig. [3] show our best fit to the 
data. Note that although the stability is determined only 
by r e /r t h and io in the two extreme limits (not visible in 
the figure), in the intermediate region of interest here 
this is not the case, so several curves are shown. Each 
solid curve segment corresponds to a single L, over the 
range of Rl tested; the dotted lines continue these curves 
for a wider range of Rl- The data are grouped into 



three inductance ranges: 6-12, 15-60, and 120-600 nH, 
indicated by circles, squares, and triangles, respectively. 
We used fixed values = 0.8, r] = 6.5, which are based 
on independent measurements, and fitted r a = 1.9 ns, 
and r c = 7.7 ns to all data. Separate values of p n vo were 
fitted to data from each of the three chips, differing at 
most by a factor of ~2. These fitted values were p n vo ~ 
1 x 10 n ^/ s ; since p n ~ 10 9 fi/m, this gives ~ 100 
m/s, a reasonable value. 

A natural question to ask in light of this analysis is 
whether it suggests a method for speeding up these de- 
vices. The most obvious way would be to increase the 
heat transfer coefficient h, which increases both I ss and 
^o, moving the wire further into the unstable region, and 
allowing its speed to be increased further without latch- 
ing. However, at present it is unknown how much h can 
be increased before the DE begins to suffer. At some 
point, the photon-generated hotspot will disappear too 
quickly for the wire to respond in the desired fashion. In 
any case, experiments like those described here will be a 
useful measurement tool in future work for understanding 
the impact of changes in the material and/or substrate 
on the thermal coupling and electrothermal feedback. 
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